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Abstract 

The algorithm-development activities at USF during the second half of 1997 have 
concentrated on data collection and theoretical modeling. Six abstracts were submitted for 
presentation at the AGU conference in San Diego, California during February 9 -13, 1998. Four 
papers were submitted to JGR and Applied Optics for publication. 

Tasks Accomplished: 

1 . Two cruises were completed. 

a. Red Tide Cruise #1 

1) August 20 - September 4, 1997 

2) Florida DEP-funded boat time 

3) Transects from the west and northwest Florida Shelf 

4) Chlorophyll a (Chi), particle absorption coefficients (ap), detritus absorption 
coefficients (ad), gelbstuff absorption coefficients (ag), remote sensing 
reflectance (Rrs) were collected at station locations, while chlorophyll 
fluorescence, attenuation coefficients (c), salinity, and temperature were 
collected underway. 

b. Florida Bay Cruise 

1) October 2-4, 1997 

2) NESDIS and NRL-funded ship time 
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3) Transects inside Florida Bay 

4) Rrs, ap, ag, chi, c, salinity, and temperature were collected. NESDIS overflew 
the bay with the VIMS Beaver equipped with the Scanning Low-Frequency Radiometer 
(SLFMR), and the Satlantic SeaWiFS Airborne Surveyor (SAS). 

2. A paper entitled ‘Satellite-Sensor Calibration, Verification Using The Cloud-Shadow Method’ 
by P. Reinersman et al. was submitted to Applied Optics for publishing. 

An atmospheric-correction method which uses cloud-shaded pixels together with pixels 
in a neighboring region of similar optical properties is described for use with high-resolution(e.g. 
30m pixels) satellite or aircraft data. This cloud-shadow method uses the difference between the 
total radiance values observed at the sensor for these two regions, thus removing the nearly 
identical atmospheric radiance contributions to the two signals(e.g. path radiance and Fresnel- 
reflected skylight). What remains is largely due to solar photons backscattered from beneath the 
sea to dominate the residual signal. Normalization by the direct solar irradiance reaching the sea 
surface and correction for some second-order effects provides the remote-sensing reflectance to 
the ocean at the location of the neighbor region, providing a known “ground target” reflectance 
spectrum for use in testing the calibration of the sensor. 

A similar approach may be useful for land targets if horizontal homogeneity of scene 
reflectance exists about the shadow. Monte Carlo calculations have been used to correct for 
adjacency effects and to estimate the differences in the sky light reaching the shadowed and 


neighbor pixels. 


3. A paper entitled ‘Semi-Analytic MODIS Algorithms for Chlorophyll a and Absorption with 
Bio-Optical Domains Based on Nitrogen-Depletion Temperatures’ by K. L. Carder, F. R. Chen, 

Z. P. Lee, S. Hawes , and D. Kamykowski was submitted to the JGR MODIS special issue for 
publication. 

This work describes MODIS algorithms for chlorophyll a concentration and 
phytoplankton and gelbstoff absorption coefficients. The algorithms are based on a semi- 
analytical, bio-optical model of remote-sensing reflectance, R rs (T), where R rs (L) is defined as the 
water-leaving radiance, Lw(X,), divided by the downwelling irradiance just above the sea surface, 
Ed(X,0 + ). The R rs (X) model has two free parameters, the absorption coefficient due to 
phytoplankton at 675 nm, a^(675), and the absorption coefficient due to gelbstoff at 400 nm, 
a g (400). The R rs model has several other parameters which are fixed, or can be specified based 
on the region and season of the MODIS scene. R rs is modeled using these parameters at each of 
the visible-range MODIS wavelengths, 7,*. R rs (7,j) is derived at each pixel from the normalized 
water-leaving radiance, L wn (Xj), measured by MODIS. These R rs (^) values are put into the 
model, the model is inverted, and a<p(675) and a g (400) are computed. Chlorophyll a 
concentration is then derived simply from the 3^(675) value. For waters with high chlorophyll a 
concentrations where R rs (412) and Rrs(443)are very low and the algorithm is insensitive, an 
empirical algorithm with R rs (488) and R rs (551) is used to estimate chlorophyll a. The algorithm 
also outputs both the total absorption coefficients, a(Xj), and the phytoplankton absorption 
coefficients, a<p(Xi), at the visible MODIS wavelengths. Since absorption per unit chlorophyll 
varies by a factor of five or more for the global ocean due to accessory pigment-to-chlorophyll 
variations and pigment packaging or self-shading, to derive accurate chlorophyll a concentrations 


based upon MODIS estimates of absorption, the MODIS algorithms are parameterized for three 
different bio-optical domains: 1) unpackaged; 2) transitional; and 3) packaged. These provide a 
variation in domain type from high-light, warm conditions to lower-light, cool conditions and can 
be identified from space by comparing sea-surface temperature to nitrogen-depletion 
temperatures(NDTs) for each domain. Algorithm errors of more than 45% are reduced to errors 
of less than 30% with this approach, with the greatest effect occurring at the eastern and polar 
boundaries of the basins. 

To partition the data into one for regions where little pigment packaging is to be expected 
but high photoprotective pigments are found (e.g., high-light, non-upwelling locations in warm, 
tropical and subtropical waters), versus locations with more packaging and little photoprotective 
pigments(e.g., eastern boundary upwelling, and non-summer, high latitude sites) two numerical 
filters were used. 

The first numerical filter compares regional field data sets to the CZCS chlorophyll 
pigment algorithm (C = 1.14 fos] 171 , r 25 = R rs (443)/R rs (551)) to check for consistency with this 
classical algorithm for Case 1 waters, which was developed with largely subtropical and summer 
temperate data. 

The second numerical filter uses the ratios ri 2 (= Rts(412)/R rs (443)) and r 25 . For waters 
with unpackaged pigments, the line r ]2 = 0.95 fos] 0 16 was used to separate high-gelbstoff data 
points from the Case 1 data. Based upon the measured a g data, the gelbstoff-rich Case 2 data had 
ag(400) values typically in excess of the relationship 0.12 [chi a] 01 . 

Those data sets generally found to be consistent with the CZCS algorithm line as well as 
occurring above the line ri 2 = 0.95 fos] 0 16 for points where r 25 > 3.0, were classified as 



“unpackaged”, and these consisted of tropical and subtropical data away from eastern-basin 
boundaries. 

There are 287 data points in this ensemble data set: 134 USF data points and 37 EqPac 
equatorial Pacific points, all measured above-surface and processed using the Lee et al. (1996) 
protocols; and 126 EqPac points, all measured below-surface using the Mueller and Austin 
(1995) protocols. 

There are 326 points in an ensemble of multi-year, multi-season data sets from the 
California Current which we label as "packaged." These consist of CalCOFI (n=303) and 
Cal9704 (n=23) data. The CalCOFI R rs data were subsurface while the Cal9704 data were 
above-surface collections. The algorithm yielded RMS1 (log fractional) and RMS2 (linear) errors 
for chi a retrieval of 0. 1 1 1 and 0.268, respectively when tested with this data set. The Type II 
RMA slope and intercept was 0.999, the bias was -0.006, and the r 2 value was 0.917 

To generate an algorithm to transition between regions and periods with packaged and 
unpackaged pigments, we developed a global data set combining the “packaged”, “unpackaged”, 
and other mixed data sets from SeaBASS. This data set has 976 data points and the algorithm 
yielded RMS1 and RMS2 errors in algorithm-derived chi a of 0.176 and 0.446, respectively. 

While algorithms appropriate for regions with packaged or unpackaged pigments can 
reduce the uncertainty in chlorophyll-a concentration from perhaps 45-50% to less than 30%, 
methods based upon space-derived data to determine when and where to apply the appropriate 
parameterization are still under development. A numerical filter has already been discussed, but 
it is only definitive for waters where r 25 > 3.0 and when no uncertainty in atmospheric correction 
exists. Also, stations with high gelbstoff concentrations can cause confusion with this method. 
For offshore oligotrophic to mesotrophic waters, however, it is a very useful diagnostic tool. 



A second space-based approach uses the fact that unpackaged pigments are usually found 
in high-light, nutrient-poor waters where small-diameter phytoplankton cells predominate. Since 
dissolved nutrients cannot be detected from space, a nutrient surrogate was sought. 

Kamykowski (1987) developed a model that explained much of the global covariance 
observed between upper-layer temperatures and nitrate concentrations. Kamykowski has since 
developed nitrate-depletion temperatures (NDTs) for the north Atlantic Ocean. These NDTs 
provide a means to observe from space a variable that indicates when and where nitrate may be 
limiting phytoplankton growth, and where upper-layer production is dependent upon recycled 
nitrogen. Such phytoplankton are typically small with unpackaged pigments and with high 
photoprotective-chlorophyll pigment ratios. 

To delimit regions of the north Atlantic Ocean with unpackaged pigments, we have 
compared sea-surface temperatures to Kamykowski’s NDTs. Figure 1 shows annual trends in 
sea-surface temperature (SST), CZCS pigment, and NDTs for the Gulf of Maine, Bermuda, and 
Barbados. The temperatures and pigments are four-year (1982-85) monthly averages from the 
AVHRR and CZCS sensors. 

Clearly, the Gulf of Maine is a lower-light, higher-nutrient environment than are Bermuda 
and Barbados, so the degree of packaging there is likely to be much higher, and the chlorophyll- 
specific absorption coefficient much smaller. By analyzing bio-optical data in the SeaWiFS 
SeaBASS archive, some preliminary functional relationships between the NDTs and pigment- 
packaging classifications for the north Atlantic Ocean were empirically derived using sea-surface 
temperature (SST) derived from the AVHRR satellite sensor: 


DOMAIN: 



1 . Unpackaged: 


SST > NDT + 3.0° C 


2. Transitional or global: NDT +1 .8° C < SST < NDT + 3.0° C 

3. Packaged: SST < NDT + 1 .8° C 

These domains for the months of May and August are shown in Figure 1 1 from the 
manuscript, based upon climatological sea-surface temperatures. 

Using AVHRR SST data from the physical-oceanographic data archive, bio-optical data 
were sorted into domains using NDTs. Data for the transition period from spring to summer 
from the NASA SeaBASS archive for the cruises MLML2, AMT4, GOMEX1, GOMEX2, and 
the North Sea were sorted into the three bio-optical domains, and the appropriate algorithm 
parameterization was applied to derive chlorophyll a values. Atlantic Meridional Transect (AMT 
4) data along 20° W longitude collected in May, North Sea data and MLML2 data collected in 
July, and GOMEX1 and GOMEX2 data collected in April and June provide a diverse set of north 
Atlantic observations that were sorted by the NDT filter and processed. The results are compared 
to those obtained by simple use of the global (transitional) algorithm. The RMS1 and RMS2 
errors for this diverse data set were 0. 1 53 and 38%, respectively, for domain-sorted data, while 
the errors grew to 0. 1 86 and 50%, respectively, when processed using global or transitional 
parameterization for the algorithm without sorting by domain. 

In summary a semi-analytical algorithm was tested using a total of 976 global data points 
from regions where the pigments were typically unpackaged, transitional or packaged with 
appropriate algorithm parameters applied for each data type. 



The semi-analytical algorithm performed superbly on each of the data sets after 
classification, resulting in RMS 1 errors of 0. 102 and 0. 1 1 1 (e.g. 0. 10 log unit), for the 
unpackaged and packaged data classes, respectively, with little bias and with slopes near 1.0. 
RMS2 errors for the algorithms were 24% and 28%, respectively. 

For the difficult transition period between spring and summer, a data set was tested that 
included the eastern-boundary upwelling region of the north Atlantic. The nitrogen-depletion 
temperature was used with AVHRR-derived sea-surface temperature to sort stations into 
“packaged”, “unpackaged”, and transitional domains. RMS2 errors dropped from 50% to 38% as 
a result of this data-sorting exercise. Since large regions of the subtropical Atlantic and Pacific 
oceans remain in the “unpackaged” bio-optical domain during all seasons and provide rather 
stable data with accuracies from 24% to 28%, it seems reasonable to expect that use of an NDT- 
based sorting algorithm with MODIS sea-surface temperatures to separate data into appropriate 
bio-optical domains will result in accuracies for the MODIS semi-analytical chlorophyll a 
algorithm that are significantly lower than our target value of 35%. 

This study is being expanded to include all oceans and a model parameterization for 
hyperpackaged pigments found in low-light, high-latitude environments has been developed for 
addition to the algorithm. 

4. A paper entitled ‘Empirical Ocean Color Algorithms for Absorption Coefficients of Optically 
Deep Waters’ by Z.P. Lee, K.L. Carder, R.G. Steward, T.G. Peacock, C.O. Davis and J.S. Patch 
has been submitted to JGR for publishing. 

Since the late 1970’s, many empirical algorithms have been developed for chlorophyll-a 
or pigment concentrations for waters from open ocean to coastal environments, but only a few 
algorithms have been developed for optical properties of the water. Applying a specific optical 



property (e.g., specific absorption coefficient), those concentrations can however be converted to 
absorption and/or attenuation coefficients. As demonstrated, in-water optical properties can be 
empirically derived from ratios of water-leaving radiance or ratios of remote-sensing reflectance 
in one step. To extend our early studies, total absorption coefficients at 440 nm (a, (440), see 
Notation for symbols used in this text) for the ocean, and those for surface pigment and for 
gelbstoff are empirically related to the ratios of remote-sensing reflectance at the SeaWiFS 
bands, which provide alternatives to the pigment-concentration algorithms for deriving in-water 
optical information. 

The total absorption coefficient of surface water, which dominates the variance of both 
remote-sensing reflectance and the diffuse attenuation coefficients, is important to many aspects 
of oceanography (e.g., water-type classification, subsurface light intensity, heat flux, etc.). 
Empirical algorithms, based on in-water measurements, have been developed for the derivation 
of the diffuse attenuation coefficient, which closely related to the total absorption coefficient. As 
the relationship between remote-sensing reflectance and absorption only weakly depends on the 
solar elevation, the spectral ratios of remote-sensing reflectance data should then be almost 
independent of solar and/or view angles. This suggests that inherent optical properties such as the 
total absorption coefficients might be derived from the ratios of remote-sensing reflectance. 
Based on smaller data sets than the present one, Carder et al. [1992] and Lee [1994] proposed an 
a,(490) algorithm using the ratio of remote-sensing reflectance at 520 and 560 nm. Since pigment 
absorption peaks near 440nm, the values of a, (440) are instead empirically related to various 
spectral ratios of remote-sensing reflectance in this study. Actually, field data show that optical 
properties at 440 and 490nm (e.g., ^(440) and A^(490), c(440) and c(490)) are all highly 


correlated. 



The pigment absorption coefficient is important for the calculation of light harvesting by 
phytoplankton pigments for use in primary production models and for determination of the 
chlorophyll-a/pigment concentration. Traditional remote sensing methods have estimated 
pigment absorption using two steps: 1) derive chlorophyll-a or pigment concentration from 
empirical remote-sensing algorithms; and 2) convert the concentration values to absorption 
values using known or assumed chlorophyll-specific absorption coefficients. That perhaps 
accounts for the different empirical and semi-analytical algorithms developed for remote sensing 
of pigment concentrations. We may, however, get the pigment absorption coefficients directly 
from the remotely measured data as briefly indicated in Lee et al. [1996a], This type of algorithm 
may be applied to a wider range of oceanic environments than the two-step algorithms, since the 
chlorophyll specific-absorption coefficient is not involved, which varies widely from place to 
place. 

Gelbstoff absorption coefficients can be used as a tracer of land runoff for coastal 
environments. No simple remote-sensing algorithm yet exists for gelbstoff absorption coefficient, 
but more complex analytical and semi-analytical approaches have been used. 

An analytical method has been developed for the derivation of oceanic absorption 
coefficients from above-surface remote-sensing reflectance. The method uses an optimization 
approach with hyperspectral data, which would take a rather long calculation time for processing 
of a satellite image. It is, however, more accurate in decomposing the total absorption 
coefficients into those of pigment and gelbstoff. For image processing purposes, more rapid 
algorithms are needed, though they may be less accurate. Recently, Hoge and Lyon [1996] 
suggested a matrix inversion technique. The technique requires, however, pre-knowledge of the 



spectral models of the inherent optical properties, and the method has not been tested with field 
data. 

In this study, empirical algorithms are developed for quickly estimating the coefficients at 
440nm for total absorption, pigment absorption, and gelbstoff absorption by directly relating 
them to the ratios of remote-sensing reflectance at two or three of the SeaWiFS wavelength 
bands. We adopt a formula of cubic polynomials using two spectral ratios. With derived 
absorption coefficients at 440 nm, the absorption spectrum of the total, the pigment, and the 
gelbstoff in the visible domain can be constructed using the approaches of Austin and Petzold 
[1986], Lee [1994], Bricaud et al. [1981], and/or Carder et al. [1991]. 

5. A paper entitled ‘Hyperspectral Remote Sensing for Shallow Waters: 1 . A Semi-analytical 
Model’ by Z.P. Lee, K.L. Carder, C. D. Mobley, R.G. Steward, and J.S. Patch has been submitted 
to Applied Optics for publishing. 

Single or quasi-single scattering theory and numerical simulations show that subsurface 
upwelling signals (SUS) can generally be expressed as a sum of contributions from the water 
column and from the bottom: 


SUS_ = SUS dp [l-e 2KH ]+ SUS B e 1KH 


(1) 


Here SUS. can be designated either the subsurface upwelling radiance (L„.), the subsurface 
irradiance reflectance (/?.), or the subsurface remote-sensing reflectance (r„); we use 
“subsurface” to mean “just below the air-water surface.” SUS dp is for optically deep waters, 



while SUS B is for the upwelling signals just above the bottom. H is the depth of the water.T K is 
usually described as an “effective attenuation” coefficient. 

The first term on the right side of Eq. 1 is the water column contribution, though in some 
studies this term was replaced by a nearby deep-water signal. The second term is the bottom 
contribution. 


Because “effective” attenuation coefficient is ambiguous and hard to determine from 
remotely measured data, and because the diffuse-attenuation coefficient for downwelling light is 
not necessarily equal to the diffuse attenuation coefficients for upwelling light, a more general 
expression for subsurface remote-sensing reflectance, r„, might be T 
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(2) 


where Kj is the vertically-averaged-diffuse-attenuation coefficient for downwelling irradiance, 
K u c is the vertically-averaged-diffuse-attenuation coefficient for upwelling radiance from water 
column scattering, and K U B is the vertically-averaged-diffuse-attenuation coefficient for 
upwelling radiance from bottom reflectance, p is the bottom irradiance reflectance, and the 
bottom is assumed to be a Lambertian reflector. 

If we define 


K= a +bb. 


(3) 



then, to first order, the diffuse attenuation coefficient K is related to /fthrough a light distribution 


function ( D ): 


K = D K. 


(4) 


In particular, 


K d = D d K, K u c = D u c K K b = D b k, 


(5) 


and we can rewrite Eq. 2 as 
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( 6 ) 


Note that ATis an inherent optical property. Thus all of the variability in the K's associated with 
the directional structure of the light field is now transferred to the light distribution functions ( D ). 

In many previous remote-sensing applications it was assumed that D d = D u c =D U B = a 
constant, which is equivalent to using a constant “effective” attenuation coefficient ( K ). 

However, studies of the sub-surface light field suggest that in general D d * D u c ?D U B * constant, 
and that the D ' s vary with the inherent optical properties of the water column. Incorrect 
assumptions for the values of the D’s will cause errors in water-depth derivations. To improve 
remote sensing of shallow-water bathymetry, for instance, how D varies with in-water optical 
properties needs to be more accurately spelled out. 



In this study, the Hydrolight numerical model was used to compute above-surface remote- 
sensing reflectance ( R rs ) values for a series of water, bottom type, sky, and depth parameters. 

From these calculated R rs values, parameterizations of D were empirically derived, and a semi- 
analytical model (SA-model) was developed for R rs of deep and shallow waters. This model is 
invertable for calculation of depths and the optical properties of the water column, where as 
Hydrolight is strictly a numerical forward model and not explicitly invertable. 

6. An abstract entitled “Methods for Utilizing Hyperspectral In-situ Light Profiles in the Presence of Wave 
Focusing and the Absence of above-water Measurements” by D K Costello, K L Carder and J S Patch is 
submitted to the Feb., 98 meeting of the AGU in San Diego for presentation. 

In simple models (e.g. black sky, smooth surface, homogenous deep water column), 
underwater light fields decay exponentially with depth. However, the complications which often 
arise in real measurements are demonstrated using hyperspectral (512 channel) light profiles 
obtained in the Florida Current. These data show amplitude perturbations in both upwelling 
radiance (Lu) and downwelling irradiance (E d ) spectra. These perturbations in L u and E d are 
sometimes temporally coherent (e.g. clouds) but are most often incoherent (e.g. wave focusing 
effects) with the greatest oscillations in E d . Furthermore, it is particularly apparent in derived E d 
data that wave focusing has not only significant amplitude but also significant spectral effects. 

2nd-order polynomial fits (P-fits) of log(Lu) and log(E d ) as functions of depth were 
generated to smooth Lu and E d amplitude fluctuations while allowing for IOP changes and 
inelastic scattering effects. The zero-order, P-fit coefficient extrapolates the spectra to the water 



surface, the first-order coefficient is equivalent to the average, spectral diffuse attenuation 
coefficient, and the second-order coefficient describes the change in attenuation with depth. 

The K u spectra (derived from P-fit-smoothed Lu data) increased with depth at the blue end 
of the spectrum and decreased at the red end with a distinct “hinge point” at 544 nm. The blue 
change is attributed to a decreasing particle reflectance (supported by ancillary c-meter data) 
while the red change is attributed to the increasing importance of water-Raman scattering with 
depth. The smoothed Ej spectra were significantly improved. Raw curves were evaluated by 
modelling Ed (z,X) to depth requiring different degrees of focusing on the direct and diffuse 
irradiance components. Note: K u and Kd spectra are used to extrapolate underwater 
measurements of Ed and Lu to the surface for comparison to satellite measurements of normalized 
water-leaving radiance. 

7. An abstract entitled ‘Characterization of Bottom Albedo Using Landsat TM Imagery’ by 
Renadette, L.A., K.L. Carder, D.K. Costello, W. Hou, D.C. English is submitted to AGU San 
Diego for presentation. 

Using Landsat TM data to characterize bottom types along the coast of the Florida Keys. 
Traditional methods of bottom interpretation, such as spectral unmixing, are not well suited to 
the limited number of spectral channels provided by Landsat. However, by applying a 
previously-developed histogram analysis method, we characterize Landsat pixels as sand, 
seagrass, or fractions thereof. This method was developed using spatial variations in brightness 
of radiance data, collected over regions of variable bottom composition. In situ measurements of 
water optical properties and bottom reflectance end members (pure bright sand, thick dark grass) 



allow us to verify our results. Our results suggest that histogram analysis of satellite data can be 
a simple and effective means of estimating and removing path radiance from a scene, and 
characterizing bottom albedos. 

8. An abstract entitled ‘Characterizations and Techniques needed for the use of Spectral 
Radiometers to Collect Plane Irradiance Measurements’ by D. C. English, R. G. Steward, D. K. 
Costello, K. L. Carder is submitted to AGU San Diego for presentation. 

Spectral plane irradiance measurements, especially downwelling plane irradiance 
measurements, are important to understanding the distribution of light in water and interpreting 
reflectance measurements both above and below the water's surface. While attempting to 
characterize radiant and irradiant spectrometers developed at the University of South Florida 
(USF) we have confronted several problems in collecting accurate spectral irradiance 
measurements, including a commercial spectral radiometer with previously understated angular 
response errors. We compare the results of several different techniques for obtaining downward 
irradiance values and present the directional response of the cosine collectors used with the 
USF spectrometers. We offer recommendations for sensor characterization and the collection of 
spectral irradiance measurements. 

9. An abstract entitled ‘Effects of Colored Dissolved Organic Matter (CDOM) and Pigment 
Packaging on Remote Sensing Reflectance Algorithms in the Southeastern Bering Sea’ by Patch, 
J.S., Carder, K.L., Lee, Z.P., Steward, R.G is submitted to AGU San Diego for presentation. 



Current remote sensing algorithms for obtaining pigment concentrations from SeaWiFS 
and other ocean-viewing sensors have been developed primarily for Case I waters. In such 
waters, the absorption due to colored dissolved organic matter (CDOM) covaries with pigment 
concentration. The southeastern Bering Sea is a region where this covariation no longer holds 
true. Regional algorithms that recognize this lack of covariance and that also take into account 
pigment packaging are required in low-light, high latitude areas such as the Bering Sea in order 
to more accurately estimate pigment concentration compared to high-light, subtropical 
algorithms. In this study, remote-sensing reflectance ratios for SeaWiFS wavelength bands are 
related to each other and to inherent optical properties (aph, ag, and fluorometrically measured 
chlorophyll a concentration) of surface waters. Shipboard data were collected during April 1996 
in the southeastern Bering Sea and the Alaskan Current as part of a NOAA sponsored Fisheries 
Oceanography Coordinated Investigations (FOCI) program. A numerical filter approach for 
identifying high CDOM and highly packaged stations based solely on Rrs ratios is examined. An 
empirical model relating ag(440) and reflectance ratios involving wavelength bands centered at 
412, 443, 490, and 555nm is developed. 

Publications: 

1. Bissett, W.P., J.S. Patch, K.L. Carder, and Z.P. Lee, 1997 “Pigment Packaging and 
Chlorophyll a-Specific Absorption in High-Light Oceanic Waters,” Limnol 
Oceanogr., 42(5), pp. 961-968. 

Anticipated Activities: 

1. The relationships between temperature anomalies and nutrients in regard to the 



packaging effect will be explored in order to reduce uncertainty in the chlorophyll algorithm 
using Bering Sea data and upwelling data from the Arabian Sea , Monterey Bay, Southern 
California Bight, and the East China Sea. 

2. Research expeditions to be completed: 

b. Lake Okeechobee 

1) January 13, 1998 

2) SIMBIOS funded 

3) Transects in Lake Okeechobee 

4) Determining stray light effects on SeaWiFS and a “dark target” evoluation of 
aerosol radiance at 412 nm. 

3. Florida Bay cruise 

1) January, 1998 

2) NOAA funded ship time 

3) Transects inside Florida Bay, SLFMR overflight. 

4. Calibration: 

1) March, 1998 

2) At U. of Arizona and Mt. Lemmon 

3) Calibrating the handheld Spectrix spectrometers. Microtops H, Reagan 
radiometers, and the Spectralon targets in the field and in the lab. 

5. Tongue of the Ocean cruise 

1) April 1-7, 1998 

2) SIMBIOS and ONR funded ship time 

3) Transects along shallow banks of Tongue of the ocean in Bahamas 



pigment (mg/m 3 ) pigment (mg/m 


4) Assessing atmospheric adjacency effects and stray light of bright/dark targets 


on SeaWiFS data. 


sst=Triangle, NDT=Diamond, chl=Star 







